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Based on specific properties of the partition function and of the quantum correlators we derive the 
exact form of the beta function in the background gauge field method for QCD with an arbitrary 
number of flavors. The all order beta function we obtain through this method has only the first two 
orders coefficients different than zero and thus is equivalent to the’t Hooft scheme. 
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I. INTRODUCTION 

The running of the gauge coupling constant with the scale is generally computed in quantum field theories using 
a perturbative approach which consists in expansion in a small parameter. In the dimensional r^ularization scheme 
beta function for QED is known at the fourth order whereas that for QCD at the fifth one [11-H- It is known that 
the first two order coefficients are renormalization scheme independent whereas the next ones depend on the specific 
renormalization procedure. In Q ’t Hooft introduced a procedure in which the beta function stops at two loops. 

In [I^we used a semi perturbative technique to compute the beta function for QED to obtain that the beta function 
stops at two loops. This method was further simplified and improved in [TTl where the exact form of the beta function 
for the Yang Mills theory has been determined from properties of the partition function and various correlators and 
by the use of the LSZ theorem. In the present work we extend our method to the more intricate case of the beta 
function of QCD with an arbitrary number of flavors. Whereas this case require more work the principles settled in 
m remain unaltered. We determine that as in [Til Ih® orders beta function stops at the two first orders 

coefficients. Note that this result is obtained without using any Feynman diagram or expansion in a small parameter. 


II. PARTITION FUNCTION FOR AN SU(N) THEORY WITH Nf FLAVORS 


We start with the gauge fixed Lagrangian for an SU{N) gauge theory with Nf fermions in the fundamental 
representation: 


^ - 9r^^d>^Al)c^ + - my)vl/, 

/ 


where, 


and, 






V 1 


Df, = 


( 1 ) 


( 2 ) 


( 3 ) 


Here is the generator of SU{N) in the fundamental representation and for simplicity we shall consider m/ = 0. We 
shall work in the Feynman gauge = 1). One can express all the fields in the Fourier space: 


Al,{x) = — ^exp[-t/c„a;]A“(A:„) 

n 

^'(x) = ^y^^exp[-ikmx]'I>{km) 

TTl 

c^{x) = 2 exp[-ikpx\c^{kp). 


( 4 ) 
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Then the Lagrangian takes the form: 

f d\C = -ii ^ klA‘^''{K)AZ{-K) + 1 ^ kle{K)c^{-K) + 

n n 

+ ^9Y.’^nA:ik„)r^^Alikm)A^''i-k^ - km) - 

n,m 

E A^^{k^)A‘^%km)Al[k,)Al[-kr.-km-kp)- 

n^m^p 

--^'^k!;)c°‘{kn)gr'’‘'Al^{km)c‘'{-kn - km) + 

n,m 

^EE fikn)j^k^n'^f{kn) + EE ^fikr,)rA‘;ik^-km)t‘^'ffikm). (5) 

f n f n,m 

The zero current partition function has the form: 

i j m fl n p 

where the exponent is considered in the Fourier space. 

It is useful at this stage to settle some of the properties of Zq. It is known that Zq apart from a factor in front is 
given by the exponential of the sum of all disconnected diagrams: 

Zq = factor X exp[E ^i] (V 

i 

where 14 is a typical disconnected diagram. Since the calculation is done in the absence of external sources all 14 
diagrams are closed and contain summations over momenta (that appear in propagators or vertices) and thus do not 
depend at all on any momenta. The factor in front is a product obtained from integrating the gaussian integrals 
corresponding to the kinetic terms. The final result has thus the expression: 

Zo = const \{{rp^ - m)^^^ expE (8) 

i j p i 

where N is coming from the Yang Mills group SU{N) and the first factor corresponds to the ghosts, the second to 
the gluon fields and the third to the fermion fields. 

We shall apply the same procedure as in m to determine whole properties of the partition function and of the 
fields. First we consider the partition function in Eq. (|6|) and introduce in the integrand the quantity to obtain: 


^o=/nnnnnn dA1;^{ki)dc''{kj)dc‘^{km)d'^fi{k„)d'i!fi{kp)exp[i / d'^xC] = 

fl i j n p m 

= I Y[Y[Y[Y[Y[Y[dA‘l;^{h)dc\kj)dc‘^{km)d'^fi{kn)d-^fi{kp)^^^exp[i j d'^xC] = 

fl i j m n p v\ ) 


=/nnnnnn dAl^(k{)dc^(kj)dd^{km)d^ fi{kn)d'^ fi{kp) 

fl i j m n p 


dAiik) 


[A“(fc)exp[f / d'^xC]\ — 


-/nnnnnn dAl{h)d^{kj)dc\km)d^fi{kn)d-^fi{kp)AZ{k)—^eM^ / 

fl i j m n p i>\ ) J 


where contains separate products over flavors and colors. 

The first term on the right side of the Eq. 

/nnnnnn dA’l{ki)dc^{kj)dc’^{km)d'i!fikn)d'^{kp)Al{k)eyL-p[i J d‘^xC]A'^{k)=+oc, 

fl i j m n p 

/nnnnnn dA‘l^[ki)dc^{kj)dd^{km)d'l!f{kn)d'^(kp)Al{k) exp[i f d^xC]A'^{k)=-oo 

fl i j m n p 


(9) 


( 10 ) 
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is zero since the e term in the kinetic term will lead to an exponential that goes to zero (see for more detailed 
explanation). Here the product satisfy the constraint: ^ A^{k)). 

The second contribution leads to the result: 

fl i j m n p 

^gk^Y.r^^A:{k)Ai{p)A^’'{-k-p) - ^gj2p^c\p)f^^A:{ky{-p -k)- 

p p 

E A:ik)Al{p)A^''{q)A^f^{-p -k-q) + 

P,(l 

- k)] X exp[i f d^xC]. (11) 

P 

We apply the same procedure to the partition function but this time introduce in the integrand the quantity , 
where / is a flavor index and r is a color one. This yields: 

Zq = -i /nnnnnn dA;ih)dc\kj)dc‘^(km)d^fl{kn)d'^Mkp) X 

fl i j m n p 

1 E + k)^{p)] exp[* J d^xC]. (12) 


Here we used the fact: 

-/nnnnnn dHyfcOc;c''(%)dc‘'(fc™)d§/z(A:„)dM//i(fcp)^’^^(fc)exp[i f d^xC]\^^,^^^^ = 0 (13) 

fl i j m n p 

since the spinors fields anticommute and there is no pairing for (there is no integration over and thus 

the result is zero. 

A similar procedure applied to the ghost field c’^{k) leads to: 

fl i j m n p 

[^fc^c‘=(fc)fc^c‘=(fc) --^gEc“(p)p''/“''''c''(fc)Aj)(p-fc)]exp[i I d'^xC]. (14) 

P 

Next we apply the operator to the Eq. ([S]) to obtain: 

^^S=/nnnnnn dAl{h)dc\kj)dc‘^{km)d^fl{kn)d'i'fl{kp) X 

fli j m n p 

^ ^ ^ A“ (fc)^ (p) (-p - k) - 

P 

-^^k>^c^{k)gr’’‘^Al{p)c‘'{-p-k) + ^^{k)-f^k^^{k)]xexp[i J d'^xC], (15) 


NfN-2{N'^ 



Zo- 


where from Eq. ([5]) we calculate: 


dki^ 


(16) 
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III. RENORMALIZATION 


In this section we shall consider all the results in section II from the perspective of renormalization. Thus the 
renormalized Lagrangian is: 

f ^ klA^%K)AZ{-K) + izi ^ kle[K)c^{-K) + 

^ n n 

+ i^Z:iggJ2KA:{kn)r^^Al{krn)A^-'{-kr. - km) - 

n,m 

A^^{kr,)A^’'{km)A‘;^{kp)Al{-k^ - km - kp) - 


n,m,p 


E Kc^(k„)9r^^K(^m)c%-kn - km) + 


n,m 

1^2 ^ ^ + ^gZ'^ Y E 

f n f n^m 


where for simplicity we drop the index r from the renormalized fields. 
Then Eq. (ITT]) will become: 


(17) 


^0=/ nnnnnn'^^M(^*)^c''(%)dc''(fc^)d^/Kfcn)d'i'/i(fcp)(-*)[-^3y4i“‘^(fc)^“(-fc)+ 

fli j m n p 

j 

^39 ^gk^ Y {p)A^‘^{-k-p)-—Z[gY {k)c%-p-k)- 


Y A:{k)Alip)A‘^''iq)A^>^i-p -k-q) + 


P,Q 


g^Z^gY'^iPh'"t‘'A‘l^{k)'i’{-p-k)]xexp[i J d'^xC]. 


(18) 


Eq. (1121) will transform to, 


Zo = -i Y[Y[Y[Y[Y[Y[dA‘;:,ih)dc\kg)dc‘^ikm)d^fiikn)d'l'fiikp) X 

f I i j m n p 

+ -^Z!^gY^ik)-f'^Ppt'^gAl{-p + k)^{p)]exp[i j d'^xC]. 

P 

whereas Eq. da yields: 

Zo = -i j X{X{X{X{X{X{dAl{ki)dc\kg)dc%km)d^ fi{kn)d-^ fi{kp) X 
fl i j m n p 

^^Zik'^^{k)k'^c^{k) - r^gZ[ Y c°'ip)p^r'’‘'c‘'ik)Al{p - k)] exp[i ^ d'^xC]. 


(19) 


( 20 ) 


Finally Eq. (TTSl) will lead to: 


^'^S=/nnnnnn dAl {h)dc\kg)dc^ {km)d^ fl (fc„ )d^ fl {k 

fl i j m n p 

i 

i72 


,[-^Ziik^A‘^^k)A:{-k) + Z,^e^{k)c%-k) + * A“(fc)rtA^(p)A-(-p- fc) - 


V '■ 

p 

-:^Zi Y kf^c-{k)gr^‘^A^p{p)c%-p -k) + fz2^{k)rk^'i>{k)] X exp[i f d^xC]. 


( 21 ) 
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IV. RELATIONS AMONG THE RENORMALIZATION CONSTANTS 


In the path integral formalism the two point gluon function has the expression: 

{n\T[Al{x,)Al{x2m) = 


lim 

T—)-oo(l—ie) 


f dApdc‘^dc‘^d'i'd'i'A^(xi)A^(x 2 ) exp[i f d'^xC] 
f dA^dc'^dc’^d’^d'^ exp[i f d'^xC] 


( 22 ) 


We apply the LSZ reduction formula in the path integral formalism and in the Fourier space: 

(0|r[A“(pi)...A^(p^)A^(A:i)...A^(fc„)]|0)^ 

/ m 

) polarizatioiifactor x const x {pi...'Pra\S\ki...kn) I 



(23) 


In m we illustrate in detail how we apply this formula to the gauge and ghost terms in the relations in Eqs. m, 
(HU), ® and (Ell)- 

LSZ formula is more intricate and complicated for fermions as it can be seen from the following equation for a 
process with two initial and two final fermions [T^ : 

out{f\i)in = {f\S\i) ~ J d'^xid''^X 2 d'^yid''^y 2 exp[-ikiyi][us'^{ki){-i^f^dfj,y^ + to)]/3i x 
exp[-ifc2?/2][M4(A:2)(-i7'"9py2 +x 

{0|T^'^2(?/2)^'/3i(yi)^ai(a:i)^a2(a;2)|0) X 


fj.xi +m)usi(pi)]ai]exp[*pia:i] X 


flX2 H- Tn^Us2 (^^ 2 )] OL2 ]exp[*p2a^2]- 


(24) 


Here ai, a 2 , /3i and (32 are spinor indices and all momenta are on shell. 

This formula is too intricate to be easily applicable to our calculations. If aj(p and 6|(p) are the operators that 
create a one particle state with charge I respectively -I one can write: 


al{p)in -t i J d^x'^{x){ij'^ d fj, + m)us{p) exp[ipa:] 

asil^out ^ i J d'^xexp[-ipx]us{p)){-n^dfj, + m)'if{x) 

i / d'^xexp[ipx]vs{p)){-i^^d^ + m)'i>{x) 


bs{p)out^i I d'^x'i>{x){i-f^^fj,+m)vs{p)exp[-ipx]. 


(25) 


Let us rewrite the first equation in (1251) in the Fourier: 

d^k 


alip) =i d'^x 


drk 


^(A:) exp[zA;x](i 7 ^^^ + m)us{^ exp[ipx] = 


(27r)" 


J J ( 27 r )4 

i4'(p)(7'"pp + rn)us(p) 


exp[ikx\'^(k){—^^k^ + m)us{p) exp[ipx] = 


(26) 


The above formula is still useless as we need to express 4'(p) in terms of the other quantities. In order to solve that 
we consider the sum: 


^ a\{p)inUs{p) = ^ i'3>{jp)in^Pp. + m)us{p)Us{p). 


Knowing that the following formula holds, 


^ = (-7>p + to), 


(27) 


(28) 







we obtain: 
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From Eq. (l29l) we extract: 


= -W{p){p'^ - m^). 


^(P) = — 



(29) 


(30) 


As stated in section II we shall take m/ = 0 in all subsequent calculations. 

Now we shall apply all these findings to Eq. (ITSl) . The results from applying the LSZ reduction formula for the 
gauge fields and ghosts are detailed in [ll|. Here we shall consider only the fermion fields. First we divide Eq. m 
by Zq which yields: 


I = terms that do not involve fermions + 


n n n n n n (fcOdc-'’(fc,)dc‘'(fc™)d^yi(fc„)d'h/KA:p) 

fli j m n p 


k)exp[i J d'^xC]. 


( 31 ) 


But the last term in Eq. dMI) is just : 

{-^)9^Z!, ^(H|r[’F(p)7^t“A“(fc)vI/(p - fc)|H). 

P 


(32) 


Then by applying Eqs. (1^ and (1501) to Eq. (15^ one obtains for Eq. (ISTl) : 


I = terms that do not involve fermions + (—*) 9 t 7 T 7^9 „ x const 

^ p2(^p f^yj.2 

5)^ S\{{p + k),s')j)Us' (pTt). 

s,s' 


(33) 


Note that in the above equation the term in brackets actually contains the vertex function 1Z“ which is known by the 
renormalization conditions. Then one can write: 

^ AT^ — 1 

HP)yjTMKhus{P + k)^ — - — gp^ ( 34 ) 

where we used the fact that for an on shell fermion p^ = m? and also as defined in the present work the vertex function 
{Vp^ = 5t“jP^)contains already a compression between two fermion states. Also note that the factor k^p^fk+p)^ 
in the limit of on shell states leads to a constant. 

Since similar procedure (see also [Ill) applies to all the fields and interaction that appear in Eqs. (ITHt . (IT^ . (l20l) 
and (j2ip these relations will become: 

I = aZs + bZ^gp^ + cZ[g^ + dZigp'^ + gZl^g^ 

Si = 52.^2 + s^Zl^g^ 
n = r2Zi +r3Z[g‘^ 

X = yZs + zZi + qZ2 + uZsgp'^+wZ[g'^, ( 35 ) 

where we absorbed all the constants in front of the terms in the new coefficients a, b, c, d, e, x, y, z, q, u, v, w, ri, 
r 2 , rs, si, S2, S3. 


V. BETA FUNCTION AND DISCUSSION 


In general in the dimensional regularization scheme similar relations exist also in other schemes) the renormalization 
constants satisfy the Slanov Taylor identities: 


2 Zsg 2 e Zig 2 e 




sV' = 






(36) 
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where d = A — e and /i is a parameter with dimension of mass. In the background gauge field method there is a great 
simplification given by the relations: 


^2 =^2 

Z 3 = Z^g = Zig. (37) 

Then one can write the four relations in Eq. (1351) in a more compact form: 

1 = (/i + /25^ + fs9'^)Zs + fiZig^ + /5-^25^ 

1 = (ti + t 2 g^)Z 2 
1 = (ci + C 2 g^)Zi 

1 = {hi + h 2 g^)Z^ + Zi{h 3 + hjig^) + h^Z 2 (38) 


From the last three equations in Eq. (l38l) we determine: 

_ (Cjti - h^ti - fesCi) + (C2ti + Cit2 - h^t2 - h^ti - fe5C2)g^ + (C2t2 - h/it2)g'^ 
^ (/ll + ^2ff^)(ci + C25^)(tl + t25^) 

From the first three equations in Eq. (1381) we compute: 

„ _ Citi + {C2ti + Cit2 - fjtl - fhCl)g^ + ( 02^2 - - /5C2)g^ 

^ (/l+/2g^ + /35'‘)(ci+C25^)(il+i25^) 


(39) 


(40) 


We shall use this last equation as a constraint. Matching the order of the coefficients with those in Eq. (IMl) we obtain 
the condition: 


(C2t2 - hit2)f3 = 0 


(41) 


from which we deduce C 2 = /14 since none of the coefficients are allowed to be zero. 

There are some coefficients in Eq. (1381) that can be determined directly form the preceding defining equations. 
These are those associated with the terms involving the gluon, fermion or ghost two point functions. Thus from Eq. 
(I2ip and the subsequent versions of it (noting that hi, I 13 and /15 are associated with the two point function for gluon, 
ghost and fermion respectively) one can compute by simple gaussian integration: 


-4(7V2 - 1 ) 
NNf - 2(1V2 _ 
2 (iV 2 - 1 ) 
NNf - 2(iV2 - 1) 


NNf -2{N‘^ -1)' 

Furthermore Eq. leads to the following useful recurrence relation: 


(42) 


1 — /13 — /15 = /ii. (43) 

Similarly from Eqs. m and (|2ni) one can determine ci = = 1 for the same reasons. By substituting in Ea. (|3^ 

the correct values for hi, / 13 , / 15 , Ci, ti the expression for Z 3 becomes: 


where. 


^3 


_ 1 + _ 

1 + vig'^ + V2g^ + vsg^ ’ 


(44) 


Ui 

Vi 

V2 

V 3 


_ t 2 
~ 1 h~ 

= C2 + 

h 2 t 2 


h3t2 


h 2 

t 2 + ^ 

hi 


hi 
h 2 C 2 t 2 


+ 02^2 + 


I 15 C 2 

hi 


h 2 C 2 

hi 


hi 


(45) 













In general in the dimensional regularization scheme the renormalization constant can be written as: 


^3 


1 + E 





(46) 


Note that similar expression exist for any renormalization scheme by making simple substitutions for e. 
The beta function is defined as: 




( 1 ) 


= -9 




( 1 ) 


dg"^ 


(47) 


We identify Eq. (IT41) with Eq. (H51) noting that the degree of divergence (given by powers in i) of the coefficients 
Ml, Ml, V 2 , Ms is zero, one or greater than one: 


[1 + + y ( fi ) g 2 + ^_ l _ g 2 + ^ 0)^4 ^ ^ ^(0)^6 + ^ ^(0)^2 ^ !^ g 2 ^ (43) 


( 1 ) 


d2) 


.,( 1 ) 


.,( 1 ) 


„( 1 ) 


,( 1 ) 


e 

To order - this leads to: 



[1 + v[°^g^] + 


,(i) 


.,( 1 ) 


„(i) 


-r + 


-9 = 


e 


(49) 


Here, 


,<“> = 4 “’+ 4 "’ + 

( 1 ) _ ih2C2t2Y°^ 


hi 


= 


hi 


(50) 


Using Eq. (1381) and the constraint C 2 = I 14 by considering simple expansion in powers of i one obtains that 

C 2 °^ = 4°^ = ^ 2 °^ = 0- Eq. (l50)l leads also to = 0- Then from Eq. (l4^ the dependence of on the coupling 
constant emerges: 


^4^ = Pog^ + Pi9\ (51) 

where /3o and /?i are coefficients independent of . We thus determined the all order shape of the beta function of 
QCD only by using global properties of the partition function and of the various two, three or four point correlators. 

According to Eq. (1471) the beta function contains only the first two order renormalization scheme independent 
coefficients and is given by: 




^3 2N ^ " 2567r4 


(52) 
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